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© Method for synthesizing analysis model and flow analysis system. 

© The present invention is applied to the design of a mold for a semiconductor or the like and provides a 
method for synthesizing analysis models, in which analysis programs corresponding to shapes and features of 
an object to be analyzed are registered in advance in a numerical analysis system, the analysis programs are 
OQmade to correspond to partial shapes divided according to the shapes and features of the object to be analyzed 
p^and are partially modelled, whereas input and output items between adjacent portions are controlled on the basis 
of the positional relationship between the portions and input and output informations of the analysis programs 
q placed in correspondence to the portions. The invention further provides a flow analysis system for computing a 
lOflow state of material to evaluate a shape of a flowpassage, a flow control condition, material and the like, the 
W system comprising input means for inputting a flowpassage shape of a mold, properties of material and flow 
O control conditions, model synthesizing means for extracting a feature of the flowpassage shape, judging an 
adaptability with the already registered analysis model from the extracted features and synthesizing flow analysis 
^models capable of analyzing the entire flowpassage, and program execution means for removing programs 
corresponding to the synthesized model to sequentially execute them. 
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METHOD FOR SYNTHESIZING ANALYSIS MODEL AND FLOW ANALYSIS SYSTEM 



BACKGROUND OF THE INVENTION 



Field of the Invention 

The present invention relates to a method for synthesizing an analysis model adapted to automatically 
synthesize an analysis model according to a shape or the like of an object to be analyzed, and a flow 
analysis system which evaluates and designs a flowpassage of a mold, a resin characteristic and a molding 
condition from a flow state of resin. 



Description of the Prior Art 

In numerical analysis for the design of products, generally, a shape of an object to be designed is 
75 mathematically modelled. According to the latest CAE system, modelling (mesh automatic division when a 
finite element method is used) of a shape-defined product is automatically carried out, for example, as 
shown in "Mold Technique - Data Book For Mold For Plastic Injection Molding" (Vol. 2, No. 11, October 20, 
1987, issued by Nikkan Kogyo Shimbunsha), Chapter 2 (Runner and Gate), pages 16 to 19. 

Further, in a conventional system, in designing a mold for housings such as electric appliances for 
20 home-use, a fluidity of a thermoplastic resin as a material is evaluated. In designing, the- matter of 
importance is an external shape as a product, which is related to a design. 

The above-described prior art has paid not attention to the following points and could not be applied to 
design a molding die for semiconductors. 

(1) Conventional model generation methods are peculiar to analysis program using analysis procedures 
25 such as a finite element method and a difference calculus, and therefore, no consideration has been taken 

in which a plurality of analysis programs are combined to form an analysis model. An analysis model is not 
formed by combining a plurality of analysis programs in a manner such that one part of an object to be 
analyzed is worked out by a difference calculus while the other is worked out by a finite element method, or 
one part is worked out by an analysis program (1) according to the difference calculus while the other 
30 worked out by an analysis program (2) according to the difference method, but preparation of a model is 
exclusively carried out by an analyst Therefore, the prepared analysis model with individual difference 
results depending on the experience of preparing models or on the difference in knowledge about the 
analysis program. It generally takes much time to prepare a model. 

(2) Semiconductor plastic packages are made of thermoset resin instead of thermoplastic resin used for 
35 housings such as electric appliances for home-use. The thermoset resin is a material which absorbs heat 

from a mold during molding to induce hardening reaction whereby viscosity is changed, complicatedly. It is 
therefore necessary to make an accurate estimate of change in viscocity of resin during molding. 

A cavity portion of a moid for a semiconductor includes therein internal elements such as a lead frame, 
chips, Au wires, etc. which constitute a semiconductor, which are possibly greatly deformed due to the flow 

40 of resin during molding. The quality of semiconductor products is directly affected by the deformation of 
these internal elements. Accordingly, it is necessary to evaluate not only the filling state of resin but also 
the influence thereof on the internal elements. However, when the internal elements are placed on a 
flowpassage of the mold which is in the form of a flattened and narrow configuration, not only the internal 
elements are deformed, due to the flow of resin, but the state of flow of the resin is greatly changed, it is 

45 difficult to develop a program for analyzing (simulation) the change. The reason is that in order to model the 
state of flow of the resin resulting from the internal elements to reflect on the program, improvements have 
to be made while adjusting to a number of experimental data, for which long periods and many steps are 
required as mentioned in the aforementioned paragraph (1). Actually, it is impossible to repeat such a trial- 
and-error program development every product. 

so (3) As mentioned in the above-described paragraph (2), in order to analyze the deformation of the 
internal elements in the mold and the like, it is necessary to perform the stress analysis and the like from 
not only the flow state of the resin but from the viscosity, flow velocity, pressure and the like obtained from 
the flow state. However, an analysis program is prepared so that data required by the analysis to be carried 
out by said program is input in a specific form, and the output data is output in the form suitable for 
expressing the result. Therefore, in the case where a plurality of analysis programs are used to perform 
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analysis it is necessary to convert input and output data. To this end. considerable *• 
Son o? programs and input and output is required as mentioned in the above-descnbed paragraph (1). 
and it is difficult to perform such analysis. 

5 SUMMARY OF THE INVENTION 

An obiect of the present invention is to provide a method for synthesizing an analysis model in which 
an a^lvsis^odercan be formed in a combination of different analysis programs with less preparation time 
,o TSeXc^ t present invention is to provide a flow analysis system which can perform analysis 
'° and J££ ? while adjusting to the shape of flowpassage flow control ^^^^ cope with 
Another object of the present invention is to provide a flow analysis system wh.ch can flexibly cope won 
the SSSSZ T shape oi flowpassage. flow control conditions, material, etc. corresponding to the change 

„ ^ %Ta^r?bS'of the present invention is to provide a flow analysis system which can perform a 
flow analysis of materials which change in characteristic according to hysteresis. 

Another object of the present invention is to provide an apparatus wh.ch can start the low ^ analysis 
w^em described above and an apparatus which can make a high-speed and accurate numerical 
Suta«on"fthe afoTsaid „ow analysis system on the host side while can make interactive process.ng 

20 tm Z^£T£^ invention is to provide a system which makes use of the aforesaid flow 
analysis system to combine and analyze analysis related to the flow analysis. ana ivzina. 

For achieving the above-described objects, the present invention prov.des a method <^ analyzing 
analyst models which comprises registering in advance analysis programs corresponding to shapes and 

25 Sures Tan object to be analyzed in a numerical analysis system, causing the analysis programs o 
correspond to partfal shapes divided according to the shapes and features in the object to be analyzed to 
^X^^VJcooHo^ input and output items between adjacent portions on he basis oT the 
po^lrrelationship between parts and the input and output informations of the analys.s programs 

30 C0 TeCsenr.rvX rt n'fu rt her provides a flow analysis system which ca.cu.ates the flow state of materia, 
to equate Sfshape of flowpassage. flow control conditions, materials, etc.. the system compnsmg input 
means ^ foMnPuffing the shape of flowpassage of -a mold, properties of material and flow control conditions 
mode? syn hesSng means for extracting features of the flowpassage shape, judging an adaptabihty > wrth 
msoect Tthe Sdy registered analyses model from the extracted feature andsynthes.zing flow analysis 

3S mSS wh ch can ana lyze the entire flowpassage, and program execution means for taking out programs 
c?rrespondr n g to the synthesized models from an analysis library to sequentially execute them. 

Moreover he present invention further provides model reconstruction means for partially "^"""V 
alreaTsTnthesized flow analysis models relative to a partial change of input data such as the shape of 

„ f ' 0W SSr tC^"S^x P e ri menta, data registration means, c—— 
defining mean^ oT materia, characteristic formula, difference evaluation formula definmg means for drtnng 
2 evaluation formula used to numerical-evaluate the difference between the expenmental data and the 
maS fha?aSeSic formula, difference evaiuation means which uses said evaluation formula i tc ^compute 
^ difference between the experimental data and the material charactenstic formula. ^ Parameter 

45 correction amount computation means which constitutes the material characteristic formula and .n which the 
ejection amount of pLmeter peculiar to the materia, determined from the experimental data is computed 
from tT Serence computed by said difference evaluation means, which means constitute materia. 

Ch3 ^~ registration means for registrating -^sis results 

so ^^fSLi^l- for displaying previous analysis result registered insaid -^^^J"* 0 " 
means and present analysis result so that both the results can be simultaneously and easily com P ared \ 

™e ^ present invention further includes know-how registration means for registrating the relationship 
between MM?***, which is a know-how for constructing a -flow analysis model and an analysis 
ST mean! ^for reaVs^ating a group of programs capable of executing an analysis model on a computer, 
55 mans pTgrams S^SSS^^^^ ™e present -invention .furthe, include* ««* interaction 
rvnthes£ n rrn?ans for dividing the shape by interaction processing; 'defining the -amount 'of feature n an 
Z£Z^u£Z£mZ allocating an analysis mode, to a partial shape, log retaining means for retaking 
ESTSSJSS^durtno interactive processing, know-how forming means for registrating a know-how ,n the 
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know-how registration means from log information to form a know-how effective for synthesizing models, 
and know-how editing means for merging know-how already registered in the know-how registration means 
and know-how formed by the know-how forming means and avoiding an occurrence of contradiction 
between the know-hows. 

5 Furthermore, for achieving the aforesaid objects, the present invention provides a work station 

comprising a bus serving as a signal transmission path between devices connected so as to start a flow 
analysis system, a bus control device for controlling said bus, a central processing unit connected to said 
bus, a main memory connected to a disk control device and said central processing unit, a display device, 
a keyboard, and a disk connectd to said disk device. 

70 The present invention further includes, in said work station in which the flow analysis system can be 
started, a communication control device which is connected to a bus to transmit data fed from the bus to 
other devices and receive data from the other devices, a bus controi device for controlling the bus, said bus 
control device effecting communication with said work station through MODEM, a central processing unit 
connected to said bus, a disk control device, a communication control device, a main memory connected to 

is said central processing device, and a work station host computer device connected to a host computer 
comprising a disk connected to said disk device. 

Moreover, for achieving the aforesaid objects, the present invention further provides a flow analysis 
system, in which the relationship between an analysis model capable of performing a higher degree of 
analysis and the output result is registered in know-how registration means, and the programmed analysis 

20 model is registered in program registration means whereby analysis related to flow analysis can be 
combined. 



25 



BRIEF DESCRIPTION OF THE DRAWINGS 



Figs. 1 and 2 are respectively views showing a schematic fiow of a processing procedure as one 
example of a flow analysis system according to the present invention; 

Fig. 3 and Figs. 4(a) to 4(d) are views showing the hardware structure and the software structure, 
respectively, as one example of the fiow analysis system; 
30 Figs. 5 and 6 are respectively graphs showing the viscosity characteristic of thermoset resin; 

Figs. 7(a) to (f) are respectively views showing the table contents of the registration units shown in 

Fig. 4; 

Figs. 8, 9 and 10 illustrate the process of forming an analysis model in the case where a 
semiconductor plastic package is taken as an example; 
35 Fig. 1 1 illustrates an internal construction of a semiconductor plastic package; 

Fig. 12 illustrates a mold used to manufacture a plastic package; and 

Figs. 13(a) to 13(h) show an image plane transition of a flow analysis system provided with an 
analysis model synthesizing method according to the present invention. 

40 DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention will be described hereinafter with reference to Figs. 1 to 13. 
First, a flow analysis system and apparatus according to the present invention will be described. Figs. 3 
and 4 illustrate the structure in hardware and the structure in software, respectively. As shown in Fig. 3. the 

45 structure in hardware comprises a work station 1 and a host computer 2, the work station 1 and the host 
computer 2 being connected by MODEMS 20a and 20b. The work station 1 is designed so that a multi bus 
11a placed under control of a bus control device 12a has a centra! processing unit 13a, a disk device 18a 
through a disk control device 14a and a MODEM 20a through a communication control device 19a 
connected thereto, and the central processing unit 1a houses therein a main memory 15a, a display device 

so 16a and a keyboard 17. The host computer 2 is designed so that a multi bus 1Tb placed under control of a 
bus control device 12b has a central processing unit 13b, a disk device 18b through a disk control device 
14b and a MODEM 20b through a communication control device 19b connected thereto, and the central 
processing unit 13b houses therein a main memory 15b. With this arrangement, data from the keyboard 17 
on the side of the work station 1 is stored in the main memory 15a through the central processing unit 13a 

55 and at the same time displayed on the display device 16. The data on the main memory 15a is stored in 
the disk device 18a through the central processing unit 13a, the multi bus 11a and the disk control device 
f y - 14a. Data is suitably transferred in such that data is transferred to the computer 2 through the communica- 
tion control device 19a and the MODEM 20a, and data from the host computer 2 is received. Data from the 
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generation unit »J^V™™^ ,„ 9 R ^ „ omental date registration unit 271. a 

5 sistirssr- s n a — u ^ - ^ • ™ 

. iUn, - - «~ * 

TSSe' now tu the shape of an object a to be analyzed and an Item a to analyze are designated I from 
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output item OUT 1 of the program Pi is converted into the input item IN 2 of the program P2 from the input 
and output items of the programs Pi and P2 by the input and output item control unit 263, and as a result, 
conversion program Hi is found. Therefore, the procedure for executing the program Hi is determined 
during the analysis of the parts A1 and A2. 

An analysis model comprising a combination of different analysis programs is formed in the above- 
described procedure. 

The synthesizing of analysis models according to the present invention has been schematically 
described. Next, a mathematical model of an analysis program will be described. 

An analysis program registered in the program registration unit 25 (see Fig. 7(f)) shown in Fig. 4(a) will 
be described with reference to Figs. 5 and 6. 

First, the fundamental matter of a model reflected in the analysis program will be described. The flow 
analysis system to which the present invention is applied is, for example, a mold for a semiconductor 
plastic package. It is necessary to handle a thermoset resin which is a package material. 

Isothermal viscocity formula for thermoset resin are expressed as follows: 



Tjo (T) = a exp (b/T) (2) 
T 0 (T) = d exp (e/T) (3) 
c(T) = f/T - g (4) 
When t=0, t> 0 (T) (5) 
When t = T 0 (T), 17 = » (6) 

where, 77 : viscosity, tjo : initial viscocity, to: gelation time, c : coefficient of viscocity rise, T : absolute 
temperature, t : time, a, b, d, e, f and g are parameters peculiar to resin not affected by the molding 
conditions. The characteristic of formula (1) at a suitable temperature T is shown in Fig. 5. 

The resin within the mold is under the non-isothermal temperature state because it flows while receiving 
heat from the tube wall. The viscocity is estimated from the isothermal viscocity formulae below: 
From the formula (1 ), 



where, u= {77/7,0 (T)} 1/cm (8) 

r = t/T 0 (T) (9) 

The formula (7) is a curve having the characteristic in which when t = 0, u = 1 , and when r = 1 1 u = «, as 
shown in Fig. 6. As shown in Fig. 6, when the state (n, Ui), that is, r = n, increases by time At and 
temperature AT from u= Ui, the state (t2, Hz) assumes. From the formula (9), t is a function of t, T, and 
an increment Arofr in the change of state can be obtained by 




1 + T 

1 - r 



(7) 




e T 



(10) 



d exp (e/T) 



AT 



(from 
formulae (3), (9)) 



T 2 



From the formula (7) 



2 




(11) 



(1-T) 



When n is changed to 7*2, the value of u is given by 
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70 



15 



hl ♦ *i + (1 2 . Xi) " ^^ (12) 

Substituting T = T 2 and u = U 2 in the following formula (13) obtained from the formula (8), 

, = ^(DU 0 ^ (13) 

gives 

By'reSCh a procedure from r = 0 to r = 1, the change of viscooity from the initial state in the 

M ^^SmwIc formulae of the aforesaid viscocity caiouiation method, the continuous formula 
^S^^co^m^ formula and the energy conservation formula. In the case of the c.rcular 

tube fiowpassage. model formulae are given as follows: 
Continuous formula: 



Q * 2* J Vk 



(15) 



20 



Motion amount conservation formula: 

i£ . 1 _± (r n 23* ) (16) 

25 as " r a vM <j r ' 



30 



Energy conservation formula: 

^^'SJSrjoonw fo? flow analysis registered in the program registration unit 25 has been 
developed on Z basis of the model formulae as described above. While in the above descnpt.on, only the 
Passage has been mentioned, it is to be noted that preparation thereof are s.m.iarly made 
40 for the case of a flat plate fiowpassage, a diffusion fiowpassage and the like. 

Now% tie present invention is described by way of the plastic flow analysis of *J^£™£^ 
c 0m £nnrt, ie tor niastic oackaqe and design of particulars of mold flowpassages and molding condrbons, and 

fele^n of^ *PP' id on the baSiS ° f the reSUlt ° f ana ' ySiS 0btain6d therefr0 " 1, 

^ processing ^ ^ be^exp^ned with reference to Rgs. 11 and 1 2 . A second uctor 

Dlastlroackage is prepared by placing a semiconductor chip 45 on which a lead frame .ndicated at 44 n 
Ra 11 is mounted nto a cavity 56 of a mold indicated at 52 in F,g. 12, and pouring a resin 53 as a plastic 
S. LterTby means of a plunger 51 from a gate 54 (an inlet) of the mold cavity 56 to harden the res.n^ 
2Ut£t£ fTow analysis to s'ee how the viscous resin 53 flows into the mold cavrty 56 ,s earned oul ; to 
« mefeby Vende" poS sibie the evaluation of particulars of mold for an object to be deseed. However, the 
Tesin material hi thermoset properties in which the resin absorbs heat from wall surfaces of me moW and 
the viscosity thereof considerably changes, which characteristic cannot be expressed by a defm te 
*eorS formufa Experimentally, a program of a flow analysis mode, with respect to a simple shape of a 
rowpasSa^ oTa mS capable oi measuring the viscocity has been developed. but the actual - construction 
« • o?a mold cavt^ Z compScated. The flow analysis mode, with such complication taken ,nto cons.derat.on 
cannoT be by o^ly one kind of calculation method. In view of this, it is necessary to Prepare an 

aSs mo'de? with complication of a fiowpassage shape taken into consider*,™ iby a common of a 
plurality of analysis programs (models). In the case where an analyst forms a model .n an attempt of 
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carrying out flow analysis with a mold cavity, registered in advance into the program registration unit 25 
stored in the disk device 18 prior to analysis are resin flow analysis programs Pi, P2 and Ps, model 
preparation procedure programs Mi, M2 and Ms peculiar to the resin flow analysis programs Pi, P2 and P3, 
respectively, and input and outpout data conversion programs H(12), H(23), H(31), H(21), H(32), and H(13) 

5 (H (i, j): a J program for converting output data of resin flow analysis program Pi into input data of resin 
analysis program Pj). A description will be made of the contents of the resin flow analysis programs Pi to 
P3 and the model preparation procedure programs Mi to M3 with reference to Fig. 8. The resin flow 
analysis program Pi is to analyze a streamlined flow which evenly receive heat termed as a circular tube 
flow from 360* around and can be applied to a flow portion of which heat source is present in a 

10 circumferential direction. Since the program Pi uses a difference calculus, a model prepared by the 
modelling procedure program Mi is-prepared by converting a shape as an object (a shape of a portion 
through which resin flows) into a circular tube having the same flow sectional area. The resin flow analysis 
program P2 receives heat only from a vertical direction termed as a flat plate flow to analyze a flow under 
the condition capable of ignoring heat transfer from a lateral direction (a low-level flat plate flow), and is 

75 prepared by use of a finite element method. Therefore, a streamline direction (a direction of flow) is 
automatically analyzed unlike the difference calculus, and need not be particularly designated and the 
modelling procedure program Ma with respect thereto is prepared by dividing a model of a flow part shape 
into triangular meshes. The resin flow analysis program P3 may be used under the condition wherein a flow 
radially proceeds from a blow opening at one point, and the modelling procedure program M3 is in 

20 accordance with the difference calculus. 

In the shape feature registration unit 242 and the shape division rule registration unit 243 of the know- 
how registration unit 24 are registered the division method of the cavity shape corresponding to the item 
(and its calculation formula) and features for calculating features of the shape of the mold cavity 41 
schematic ally shown in Fig. 9. Here, the item for calculating the feature of the cavity shape and the content 

25 of the shape division method wili be described. Main factors which govern a hardening reaction of resin 
include a heat transfer action from wall surface of a mold and internal elements (such as chips). It is 
therefore necessary to divide a shape at a point at which the heat transfer state in the mold cavity is 
different, that is, at a point in which the resin flow is separated by the internal elements 42 and then flows to 
obtain a portion in which the heat transfer state is constant. A criterion for the division of shape is a shape 

30 and a position of the chip 42 forming an internal element. Considering now the ratio between the width g of 
cavity 41 and the width h of chip 42, in the case 'where the width h of chip 42 is substantially the same as 
the width- g of cavity 41 , it is regarded that the flow is branched into upper and lower flows; conversely in 
the case where the width h of chip 42 is smaller by a certain ratio than the width g of cavity 41, it is 
considered that the flow i$ branched into four upper, lower, left and right flows. If the ratio a ( = h/g) between 

35 the width h of the chip 42 and the width g of the cavity 41 is registered in the shape feature registration unit 
242 and the description of judgement whether the division is divided into two upper and lower portions of 
the chip 42 by the value of said ratio a or into four upper, lower, left and right portions of the chip 42 is 
registered in the form of IF to THEN in the shape division rule registration portion 243, the division of the 
shape of the cavity 41 according to the shape feature becomes possible. In addition to the above, the shape 

40 feature forming the criterion of division includes the ratio between the length b of chip 42 and the length 1 
of the cavity 41. 

While the feature value forming a criterion for allocating analysis programs on a divided portion in 
addition to the feature item forming a criterion for division of the shape of cavity is registered in the shape 
characteristic registration unit 242, items for deciding the analysis programs Pi, P2 and P3 include a 
45 streamline direction (a factor for deciging programs Pi and P3) and a shape (a height of flow forming a 
factor for deciding program P2, uniformity of a heat source for deciding program Pi, that is, a 
longitudinal/lateral ratio of shape). Among them, the feature forming a criterion for judging a streamline 
direction considered is a ratio (g/i) between a width g of cavity 41 and a width (diameter) i of gate 43. 
Where the ratio g/i between the width g of cavity 41 and the width i of the gate 43 is large, it can be 
50 regarded to be a flow from a narrow place to a wide place, and the streamline can be judged to be radial. 
Conversely, where the ratio g/i is small, the streamline can be judged to be parallel. In this manner, the 
• feature items forming a criterion for deciding an analysis program and formulae thereof are registered in the 
shape feature registration unit 241, and feature judging criteria (a, 0, etc.) with respect to the feature 
calculation result forming a criterion for deciding an analysis program are registered in the form of IF to 
. 55 THEN in the shaft feature judgement rule registration unit 246. Now, the processing procedure shown in 
- ti Fig. 1 will be described. When an analyst operates the keyboard 17 to input "Flow analysis" as an analysis 
* i : t Htem together with a name "X" of mold cavity as an object to be analyzed, the inputted analysis demand is 
* taken into the system by the input unit 22 (100), and the name "X" of mold cavity is transmitted to the 
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control unit 21 whereby particulars of cavity shape with respect to the "X" wrthin the disk dev,ce 18 s 
Searched Irough the shape registration unit 241 by the control unit 21. In the case where «™ analysis 
moTercannofbVsynthesized merely by information registered in the know-how registration un. 24 the 
mode intact on synthesizing unit 267 of the model synthesizing unit 26 is started by the control unrt 21. 
T„^tLl^r^2l2rr?npul necessary information through the output unit 23 to complete a flow 
TnlsTs Z£*ZZ?£Z f informal being processed is retained by the log retainer unit 268, and 
u D S ^ t^nation processing by the mode, interaction synthesizing unit 267. the information .s converted 
Ttie Sow foling unit 269 into the form capable of being registered into the know-how registraton 
V^Z^oLJ** thereof with respect to the existing know-how is solved to register it .nto the 
, Sow-how registration unit 24. Since in this example, particulars of cavity shape of the mold wrth respect to 
*e"x« are nofregistered, the mode, interaction synthesizing unit 267 is started by the control unit 21 and 
*e input Z shape particulars is demanded through the output unit 23. on the bas.s of wh.ch the ana.yst 
SJtsSe shape particulars of the mold cavity with respect to the "X". Then the shape P^">" 
eSned^e log retainer unit 268 and function similarly to that being registered m the shape region 
, unrt 241 That!s when the shape registration unit 241 is referred to. the log retainer unit 268 .s also referred 
to n 0) ™e^he contro. unrt 21 instructs the shape feature calculation unit 260 to emulate he 
feature of the input shape. On the basis of this instruction, the shape feature calculat.cn unrt 260 reads the 
Xaoe ftfcre item and the calculation formula from the disk device 18 through the shape feature 
egTsL^ unit 242 and further reads the shape particu.ars with respect to the "X» from the 
registration un,^ ^ calculates the feature values correspond.ng to the 

! ° Se i ems aSS fS tion result registered in the shape feature registration unit 242 (120). .n this 
manner after tS features have been calculated, the shape division unit 264 reads the feature values 
reSe edln the shape feature registration unit 242 and the shape division rules registered in the shape 
Zstn aite regUation unit 243 on the bais of instructions for dividing the shape of the mold cavity from 
* me contro unrt 2 7 so as to judge how the shape of the mold cavity is divided. As a result, the mold cavrty 
T^Te^TpoZs A to D, for example, as shown in Fig. 10 (130). On the basis of the aforesaid I resurt 
of division the shape division unit 264 calculates the shape particulars of the d.v.ded portions A to D from 
L sTape pSiculafs of the moid cavity registered in the shape registration unit 241 . and names of parts 
and date of shape particulars thereof are delivered to the partial shape reg.strat.on un.t 245 (140). Then 
30 Scut^ ol part^apes with part names indexed are registered in .the disk device 18 by the parhal 
sSS regist aSn unit 245. The shape division unit 264 further calculates the positional *» 
XSAto D, which are registered in the divided portion position registration unrt ^50)_A 
shown in fig. 10 as an example, the divided portions B and C are positioned .n parallel next tc the drjded 
portion A, and the divided portion D is positioned next to the divided portions B and C. Thus, the pos.t,onal 
reiationshio of the divided portions A to D are registered. 

Xr C ^vision of shape, the control unit 21 instructs the shape feature calculation unit 26 I to 
calcutS *e features which form a criterion for allocation of analysis programs with respec to the divided 
Mrtions Ato 0 on the basis of which the shape feature calculation portion 260 reads partial shapes of the 
dieted portions A to D registered in the partial shape registration unit 245 and the feature items and 
40 cacumT^\i Leof registered in the shape feature registration unit 242 to calculate the feature 
" "SL i^a coTdtnce with the' calculation formula, which values are then registered in 

reaistration unit 242 (160). Then, the control unit 21 instructs the partial shape feature judg.ng unrt 265 to 
StTlJSii h. tJ» of allocation of analysis programs of the divided portions A to D. The partal 
shape feature judging unit 265 reads the rule in the form of IF - THEN which presume factors for deciding 
4S he ana7y7s programs from the feature values from the shape feature judging rule registration unit ^246 and 
he feature values of the divided portions A to D from the shape feature <*^^™%™ TJojuI 
which the factors forming a criterion for allocation of analysis programs of the d.v.ded P°* ,ons ° p ™* 
presumed result is registered in the partial shape registration unit 245 (170). For example with respect to 
Sie d^ded Norton A. a stream line "radial" and a flow height "height" are judged by the rafco between the 
widtS o? gate which is an inlet and the width of flow; with respect to the divided portions B and a flow 
height "low" is judged; and with respect to the divided portion D. a stream l.ne "parallel" and a flow he.ght 

"^tfte/thVteature judgement with respect to al. the divided portions A to D has been terminated in a 
manneras described above, the control- unit 21 instructs the analysis program allocation unit 26 to allocate 
ZSysi programs ^th espect to the divided portions A to D. The anaiysis program allocation unit . 261 
%£*E££7** divided portions A to D registered in the partial shape registration unrt 245, which 
enures are registered in the program registration unit 25. Among the flow analysis programs those 
matched in program using condition to the features of the divided portions A to D and names of model 
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preparation programs peculiar thereto are registered in the partial shape registration unit 245 (1 90). As the 
result, the analysis program P 3 and the model preparation program M 3 , the analysis program P 2 and the 
model preparation program M2, and the analysis program Pi and the model preparation program Mi are 
allocated in the divided portion A, divided portions B and C, and divided portion D, respectively. Thereafter, 
5 the control unit 21 instructs the model preparation unit 262 to read the model preparation programs of the 
divided portions A to D registered in the partial shape registration unit 245 and particulars of shapes thereof 
so as to prepare models for each of the divided portions A to D (200). The thus prepared models have their 
names of parts registered as indexes in the disk device 18 through the model registration unit 240. 

After preparation of models with respect to the divided portions A to D has been terminated (210), the 

to control unit 21 instructs the input and output item control unit 263 to control the input and output items 
required in continuously analyzing the divided portions A to D. The input and output item control unit 263 
reads the positional relationship of the divided portions A to D from the divided portion positional relation 
registration unit 244, and names of programs to analyze corresponding to the divided portions A to D from 
the partial shape registration unit 245 so that data conversion programs for controlling the input and output 

15 data of a connecting section of analysis programs for which input and output item control of analysis data is 
required, that is, the divided portions A — B, A — C, B — D and C — D are searched by the program 
registration unit 25. As the result, A — B corresponds to data conversion program H32, likewise A — C to 
data conversion program H32, and B — D and C — D to data conversion program H21, respectively. These 
conversion program names are registered in the divided portion postional relation registration unit 244 

20 (200). Thereafter, the control unit 251 instructs the analysis model consolidation unit 266 to consolidate 
analysis programs and models corresponding to the divided portions A .to D, respectively, and input and 
output conversion programs so as to prepare a single model as a whole (230). In the analysis consolidation 
unit 266 are read analysis program names corresponding to the devided portions A to D from the partial 
shape registration unit 245, modelled results corresponding to the divided portions A to D are read from the 

25 model registration unit 240, and data conversion program names corresponding to adjacent portions are 
read from the divided portion positional relation registration unit 244, which are then edited as the analysis 
execution procedure. In the program execution unit 28, according to the aforesaid analysis execution 
procedure, a flow analysis model as a whole within the moid cavity is automatically formed in the form 
wherein analysis programs Pi, P2 and, Pa are mixed. 

30 As described above, according to the present embodiment in the case where analysis for which 
analysis procedure is not established as in the flow analysis within the mold cavity in connection with the 
designing of a semiconductor plastic package is carried out, a plurality of existing analysis programs are 
combined to automatically form an analysis model which enables a flow analysis as a whole. 

Next, how the analysis model synthesizing function provided in the flow analysis system according to 

35 the present invention behavior will be described with reference to the image plane transition shown in Fig. 
13. 

This example shows the process in which an analyst inputs a proposed design of a mold for a 
semiconductor to synthesize flow anaylsis models according to the shape or the like thereof. 

First, a designer instructs to carry out flow analysis in connection with particulars of a mold flowpassage 

40 registered in advance. Then, as shown in Fig. 13(a), an image plane is displayed. The image plane is 
composed of a message display window 61 for displaying a message from the system, a window 62 for 
displaying the entire figure of a mold which always displays the entire mold flowpassage, a window 63 for 
displaying a partial view of a mold which displays an enlarged view of a part of a mold flowpassage, a 
figure in which a point of view is changed or a part of a mold after divided, a menu forming a criterion when 

45 ' a user divides a shape in an interactive manner, for example, a branch type menu 64 for selecting a flow 
branch state of resin in a mold, and a shape type menu 65 for selecting a basic shape in which 
corresponding relationship with an analysis program is obvious. 

First, an analyst divides a mold shape indicated in Fig. 13(a) into units that may be analyzed by an 
analysis program, At that time, interactive division/automatic division selection areas (69, 70) defined on the 

so mold entire figure display window 62 are selected by used of a mouth or the like to select modes, whether 
the shape division is effected automatically or interactively. In the case where the automatic division is 
selected, shape features in terms of shape division registered in advance are calculated, and the shape 
division according to the aforesaid feature is automatically carried out. The moid shape is divided into units 
of the basic shape type which is apparent in the corresponding relation with the analysis program as shown 

55 in Fig. 13(c). Fig. 13 (d) shows that an enlarged function 71 defined on the message display window of Fig. 
13(c) is selected by use of a mouth or. the like, and the result of the shape division is enlarged in the 
periphery of the moid cavity., Jt is understood that the shape division corresponding to the case where 
upper, lower, left and right branch flows of chips occur in consideration of the branch flow of resin caused 
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bV tSSS? S r ^l the interactive division is selected, the user uses the branch 
type menu £ or shTpe type menu 65 as shown in Rg. 13(b) to divide the mold shape wh,le grasping the 
feature of shape in an interactive manner. In the stage in which the shape division has been terminated, an 

e imsae olane as shown in Fig. 13(c) results similar to the automatic division. ^ 
Afte TeTaoe^Msion has been terminated so that the mold shape has been divided into un.ts that 
may JSjTbJ "he anlysis program, suitable analysis programs are allocated to the >J~£*~ 
oortonS to tteanZsis program allocation step, an image plane having a structure shown in Rg. 13(c) is 
Saved T^e Plane is composed of an actual shape display window 66 for displaying an art* 

,o 22* ■ M^TcLed in shape, a window 67 of a list of analysis programs ^r d.splay.ng anaK/sis 
oroo^ narTet prepared by types of resin flows such as a circular tube flow, a flat plate flow and a 
5SSS ; 0 w^d P TnTode. sC disp.ay window 68 for displaying mode, shapes provided on the d,v,ded 

POrti F^ f a a dTv°d d ed portion (in this case, a runner portion in a mold flowpassage) as displayed on the model 
» shape dspCwSdSw 68 of Fig. 13(e) is displayed. Next, the value of the shape feature of the divided 
t^VSoMrtm^*^ is judged from the shape feature value, and an analyse program 
S^S^JTS^i is decided, In the case of Rg. 13(e). the type of the analyse program is 
decSed to be a circular tube flow from a flatness of the divided portion having a shape of a sem,c,cula 
hnSm NexY a modelling program peculiar to a circular tube flow program is started, details of a model 
20 S5T£ o e ^circular tube are decided from the particulars of the divided portion having a shape of a 
semickcu^om. an^ a mode, shape as shown in Rg. 13(f) is displayed. The above-descnbed aHocation 
of modern caSed out for all the other divided portions, and finally, analysis models wrth respect to the 
ia^SSu^SX Rg. 13(g) can be synthesized. An image plane shown in Fig. 13(h) is an , enlarged 
view of *e Periphery of a cavity among mode, shapes with respect to the entire mold shown ,n Rjl3g. 
25 case of Ws Sample, it I understood that a flow model within the cavity is modelled as a diffuse 

ITp^w Tin the vicinSy of an inlet, as a flat plate flow program at the upper, lower and rear portions of 
the chiD and as a circular tube flow program sideway of the chip. 

in this manner the analysis programs allocated to the divided portions are executed while beng 
seqlSy oonnected in reliance'upon the positional relation of the respective portions registered at the 
™ Km! nf th* <;hane division whereby the flow analysis of the entire moid can be made. 

As descled to the present embodiment, analysis models which enable the resm flow 

ana lCL of *e entire mold flowpassage including the interior of the mold cavity in connection wrth he 
d^o?.^!SdS5 plasSc parage can be simply synthesized as a combination of bas.c flow 

3S ^^An^xample of the case where the flow analysis according to the present invention is applied to the 
eval^tio! I of materia, which evaluates viscosity of resin in the flowing state withm the mold will be 

d8S RrTthTtf 2 7'for presuming parameters of resin characteristic formula which is a constituent element 

of the oresent system will be briefly described. ,, n t => 

^ A semiconductor plastic package as an object in the present system is sealed by and formed of a 
^o^T^^l^J^ reaction by heat. The basic formu.a (see formula (1,) indicative of 
thHsco ityTha'aSeristte of resin has six parameters a, b, d, e, f. and g perculiar to the resm. In prior art 
ttey hav bee decided in a trial and error manner while manually depicting a graph. Therefore ,f rt takes a 
o«Id moS f than Tone month in order to decide fully satisfactory parameters, accuracy of the v,scoc,ty 
45 ^L ^^Z,eTZe^ng on the deciding person. The quality (good or bad) of the parameters of 
ST2iSST^«ipenSi on whether or not the apparent mean viscocity obtained through simoon 
b7u"r o Hormula (1) with the decided parameters coincides with the apparent means viscocity obtained 
from a measuring device. The viscocity is obtained in the following manner. 

An TorTevaluated from an apparent mean viscocity ^calculated by use of the formula of Hagen- 

so Poiseuille below: 

TheV^ apparent mean vfscocity. AP : pressure loss. 0: shape resistance. Q : flow rate 



where no: apparent mean vi&uu*-ny, ^ • ~ — ^ — . - . . . rtH _ e tho 

^om the flow pressure loss AP of the measured data, and an apparent mean vscocity „o obtained as the 
/riutt of SuLon under the same condition as the measuring condition by use of formu.ar (1). The 
evaluation formula used is 

a viscocity formula of (a. b, d. e. f. g) in the form Rector r 1/!^^ 
diameter, temperature of mold). d( ?. ™) : distance between and ,o. ex ( , - „o) 2 . | v M • To obtain the 
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optimum parameters of resin-viscocity formula is to obtain a = ao in which G(a) is minimum. That is, 
G (or 0 ) = m ^ n G (a) (23) 

ao given by the formula (23) is to be obtained. Further, if F (a) is defined by 

FCcn - J ¥^ ) ••• C24) 

that is , (f x («) , --->f 6 (#)) T 



the formula (23) corresponds to obtain a = ao as given by 

F (a) = 0 (25) 

If F (a)is teller developed around a, 

*F (a +Aa) = F (a) + F (a)*Aa+" # (26) 

Ignoring terms higher than second order of Ao, 
Act = -[ F (a)r u Aa (27) 

The repetition formula of formula (25) of a non-linear equation can be defined as follows. 

or n = a n + Aot n 

= a„-[ F (cr n-1 )]~~ 1 * F («„.,) (28) 

By providing a suitable initial value of a n , parameters of a viscocity formula can be determined by repeating 
the repetition formula shown in the formula (28) till it is converged. 

The above-described argonism is realized by the material characteristic formula parameter assuming 
unit 27 (see Fig. 2). In the case where an analyst uses the present system to presume the viscocity formula 
parameters, the experimental data registration unit 271 of the unit 27 is started by the control unit 21 in 
accordance with the input of the analyst from the keyboard 17, and the flow pressure loss P obtained from 
the measuring device is registered in the disk device 18 through the communication circuit. When the flow 
pressure loss is registered in the disk device 13, the experimental data registration unit 271 uses the 
formula (21) to change the apparent mean viscocity, and the data is also registered in the disk device 18 
(500). When the analyst inputs the instructions from the keyboard* 17, the characteristic formula defining unit 
272 of the assuming unit 27 is started by the control unit 21. The characteristic formula defining unit 272 
demands the analyst to define the viscocity formula (formula (1)) through the display device 16. When the 
analyst inputs the viscocity formula (formula 1) in the keyboard 17 in accordance with the meassage of the 
disk device 16, that viscocity formula is taken into the system (510). Next, the control unit 21 starts the 
difference appraisal formula defining unit 273 demands the input of the difference appraisal formula (formula 
(22)) through the disk device 16. When the analyst inputs the difference appraisal formula (formula (22)) 
from the keyboard 17 in accordance with the message of the display device 16, the difference appraisal 
formula is taken into the system (520). Then, the control unit 21 starts the difference appraisal unit 274. The 
difference appraisal unit 274 uses the difference appraisal formula to calculate the difference as follows. 
First, the measuring condition stored along with the measuring data in the display device 16 is taken out, 
and sent to the model synthesizing unit 26 through the control unit 21 to form a flow analysis model 
capable of being subjected to flow analysis under the same condition as the measuring condition. Then, the 
program execution unit 28 executes the analysis program in accordance with the flow analysis model. At 
that time, As a viscocity formula, the formula inputted by the characteristic formula defining unit 272 is 
used. Using the apparent mean viscocity at various parts and time obtained by the execution of the 
program execution unit 20, the difference is calculated by the didfference appraisal formula defined by the 
difference appraisal formula defining unit 273 (560). If the difference is less than an allowable value, the 
difference appraisal unit 274 terminates. However, if the difference exceeds an allowable value, the control 
unit 21 demands the parameter correction amount calculation unit 275 to start (570). The control unit 21 
starts the parameter correction amount calculation unit 275. The parameter correction amount calculation 
unit 275 obtains a correction amount of a viscocity formula parameter by the repetition formula shown in the 
formula (28) to correct the viscocity formula parameter (580). After correction, the parameter correction 
amount calculation unit 275 demands the control unit 21 to start the difference appraisal unit 274, then the 
control unit 21 starts the difference appraisal unit 274. The difference appraisal unit 274 repeats the above- 
described process to make the difference less than the allowable value. 

In the manner as described above, the viscocity formula parameters are obtained. 
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1 20 Calculate feature 
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Fig. 2 
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Fig. 3 
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Fig, 4(a) 
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24 Know-how registration unit 

25 Program registration unit 

26 Model synthesizing unit 
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28 Program execution unit 



Fig. 4<b) 
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240 Model registration unit 

241 Shape registration unit 
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Fig. 4(c) 
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Fig. 8 
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4. In a system for calculating a flow state of material to evaluate a flowpassage shape, a flow control 
condition, material and the like, a flow analysis system comprising experimental data registration means for 
registering experimental data, characteristic formula defining means for defining a characteristic formula 
representative of material characteristic* difference appraisal formula defining means for defining an 

5 appraisal formula in order to numerically evaluate a difference between experimental data registered by said 
experimental data registration means and the value of the characteristic formula defined by said characteris- 
tic formula defining means, difference appraisal means for computing the difference between the experi- 
mental data and the value of the characteristic formula using said appraisal formula, parameter correction 
amount computation means which constitute a characteristic formula and in which the correction amount of 

70 parameters peculiar to the material decied from the experimental data is computed from the difference 
computed by said difference appraisal means, which means constitute material characteristic formula 
parameter presuming means to be added, wherein the flow analysis with respect to material whose 
characteristic changes with hysteresis is carried out. 

5. The system according to claim 2, which includes analysis result registration means for registering 
75 analysis results obtained by said system, and result display means for displaying the previously analyzed 

result registered in said analysis result registration means and the present analysis result so that both the 
results can be simultaneously and easily compared, and the properness of the analysis result can be 
evaluated. 

6. The system according to claim 2, which includes know-how registration means for registering the 
20 relationship between the shape feature which is a know-how for constructing an analysis mode! and an 

analysis model, and program registration means for registering a group of programs capable of executing 
an analysis model on a computer, and in which an object to be analyzed and a range of analysis can be 
sequentially extended. 

7. A flow analysis system comprising model interactive synthesizing means for dividing a shape by 
25 interactive process, defining a feature amount in an interactive manner and allocating an analysis model to a 

partial shape, log retainer means for retaining input information (log) during interactive processing, know- 
how registration means for registrating the relationship between a shape feature which is a know-how for 
constructing an analysis model and an analysis model, know-how forming means for forming a know-how 
registered in said know-how registration means, said know-how being effective for synthesization of models, 

30 know-how editing means for merging the know-how already registered in said know-how registration means 
and the know-how formed by said know-how forming means and avoiding occurrence of contradiction 
between the know-hows, and model synthesizing means for synthesizing flow analysis models capable of 
analyzing the entire flowpassage, wherein even with respect to an object in which flow analysis models 
cannot be synthesized in terms of the relationship (know-how) between the already registered shape feature 

35 and the analysis model in said model synthesizing means, the flow analysis models can be synthesized by 
once subjecting them to interactive processing. 

8. The system according to claim 2, comprising a work station composed of a bus serving as a signal 
transmission path between connected devices, a bus control device for controlling said bus, a central 
processing unit connected to said bus, a disk control device, a main memory connected to said central 

40 processing unit, a display device, a keyboard and a disk connected to said disk device, wherein said main 
memory stores said model synthesizing means and said program execution means, said disk stores the 
adapting relationship between the shape feature and the analysis model and an analysis library, and said 
central processing unit causes said main memory to be started. 

9. The system according to claim 8, wherein a communication control device connected to the bus in 
45 said work station to transmit data fed from said bus to other devices and receive data from the other 

devices is provided, and there is provided a work station and host computer device connected to a host 
computer comprising a bus for carrying out data communication with said work station through MODEM, a 
bus control device for controlling said bus, a central processing unit connected to said bus, a disk control 
device, a communication control device, a main memory connected to said central processing device, and a 

so disk connected to said disk device, wherein a memory on the side of said host computer stores said model 
synthesizing means and said program execution means, said disk stores an analysis library, said input 
means and said model synthesizing means are executed on the side of said work station, the synthesized 
flow analysis model is transmitted by said communication control device, said program execution means is 
executed at high speed on the side of said host computer, the result of computation is transferred by said 

55 communication control device to said work station, and the result thereof can be confirmed on the display 
on the side of said work station. i 

10. The system according to claim 6, wherein output data obtained from the result of flow analysis is 
used to register the relationship between the analysis model capable of performing analysis of higher 
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